Stroboscopic Generation of Topological Protection 
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Trapped neutral atoms offer a powerful route to robust simulation of complex quantum systems. 
We present here a stroboscopic scheme for realization of a Hamiltonian with n-body interactions on 
a set of neutral atoms trapped in an addressable optical lattice, using only 1- and 2-body physical 
operations together with a dissipative mechanism that allows thermalization to finite temperature 
or cooling to the ground state. We demonstrate this scheme with application to the toric code 
Hamiltonian, ground states of which can be used to robustly store quantum information when 
coupled to a low temperature reservoir. 

PACS numbers: 03.67.Pp, 03.65.Vf,03.67.Lx 



Introduction.- Among the most exciting aspects of quan- 
tum simulation is the possibility of generating and study- 
ing exotic quantum phases such as those possessing topo- 
logical order that can be used to robustly store and pro- 
cess quantum information. The Hamiltonians governing 
these phases frequently require more-than-2-body inter- 
actions that are hard or even impossible to realize nat- 
urally. This difficulty has spurred much theoretical and 
experimental effort in the artificial engineering of Hamil- 
tonians, particularly for trapped neutral atoms [l|. Many 
proposals have been made for the generation of 2-body 
Hamiltonians using static emulation schemes and some 
experimental realizations have appeared @, Q • Specific 
proposals have appeared for generating n-body interac- 
tions 0] , but have focused on static emulation. 

We present here an alternative, dynamic emulation ap- 
proach to systematic generation of ?i-body interactions 
that is based on sequences of control pulses which in- 
dividually realize 1- and 2-body operations on internal 
atomic levels. We show that this stroboscopic realization 
of the Hamiltonian can be implemented simultaneously 
with a dissipative thermalization protocol to stabilize the 
system from the effects of imperfect quantum operations 
and environmental noise. In the zero temperature limit, 
this can be view as replacing algorithmic error correction 
in an equivalent quantum circuit model with a dissipa- 
tive procedure to remove errors [B| . The resource require- 
ments for this thermalization protocol are different from 
those of algorithmic error correction, and may be more 
accessible to experiment in the foreseeable future. We il- 
lustrate the approach here with stroboscopic generation 
of the 4-body toric code Hamiltonian, which constitutes 
one of the simplest exactly solvable models with a ground 
state topological phase 
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tice. The ground state of this model possesses topolog- 
ical order, and therefore has anyonic quasiparticle ex- 
citations and, on a lattice with periodic boundary con- 
ditions, an emergent topological degeneracy. Quantum 
information can be encoded in this ground state degener- 
acy and manipulated with controlled creation and braid- 
ing of anyons d, 0] ■ In a finite sized system [H, Q , the 
topological order of the ground state and gap to excited 
states protects against decoherence and loss of quantum 
information due to noise provided the system is coupled 
to a low temperature bath. Our analysis below will pro- 
vide a scheme for generating both Hq and an effective 
low temperature bath, realizing the topological protec- 
tion characteristic of the toric code. 

The physical context for our analysis is a set of ~ 250 
individual 133 Cs atoms trapped at the sites of an address- 
able simple cubic optical lattice [loj ]. A lattice spacing 
of 5 (im allows essentially perfect addressability [H[ . 
The orbital degrees of freedom are frozen on the time 
scales relevant to our analysis and we need consider only 
internal atomic degrees of freedom. Two hypcrfine levels 
(e.g., \F,mp) — |4,4) , 1 3 , 3) ) define a 2-level pseudospin 
system. We realize Hq in the interaction representa- 
tion defined by the pseudospin energies. Auxiliary in- 
ternal levels are used to realize 1-spin and 2-spin quan- 
tum operations, using optical frequency Raman pulses to 
generate arbitrary single-spin operations and excitation 
of one atom to a Rydberg state, e.g., the n as 80 state, 
to generate controlled-phase gates, CPHASE [12|. To 
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where <jj denotes a Pauli operator on the links of a square 
lattice and v/p denote the vertex/plaquette of the lat- 



achieve thermalization or cooling, the Hamiltonian Hq c 
is supplemented by coupling the primary system spins 
to an ancillary set of pscudospins that will be dissipa- 
tively controlled to simulate a thermal reservoir. Since 
the pscudospins are localized at the sites of a cubic lat- 
tice, one can choose to either realize Hq C on a single 
plane using a surface code [HI or in a three-dimensional 
cubic array with toroidal boundary conditions realized 
by SWAP operations (see Fig. [T]). 

Effective Hamiltonian Evolution.- Given the ability to 
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Figure 1: (Color online) An embedding of the toric code 
into a cubic lattice. Twisted-periodic boundary conditions are 
imposed by SWAP-gate shuttling along auxilliary sites, indi- 
cated by dashed-arrows. Bold lines connect logical nearest- 
ncie hbors. 



perform both Rydberg-induced CPHASE gates between 
atoms in neighboring sites and arbitrary 1-body rota- 
tions, exp (—iOcTj), on individual atoms, where 9 is a vari- 
able phase angle, sequences of these operations can be 
chosen to generate effective n-body interactions through 
high-order terms in the Magnus expansion [lij . allowing 
stroboscopic simulation of a broad class of Hamiltonians. 
Consider the operator sequence, U n U n -\ ■ ■ ■ U2U1, where 
the Uj are the 1- or 2-body gates described above. Effec- 
tive interactions are found through: 

H cff (t) = j\n(U n U n - 1 ---U 2 U 1 ) 

= J ^ Uj - £ § [ In U h \nU h ] + o(\\hxU\f 

3 ' 3<k 



Consider now simulation of the 4-body interactions in 
Hj c . We use the notation Uj{4>) = e _i ^ Ej and define 



Si = a z a y a°a°, 



£ 2 = a°a x a v a°, and £3 = a°a°a x a z , 



where er° is the identity operator. For simplicity, it is 
assumed that each Uj(a) takes a time r to execute. We 
construct the operator sequence, 

(a, /3, 7) - U 12 (a,/3)U 3 (j)Ul 2 (a,/3)Ul(j), (2) 

where U 12 (a,P) = U 2 (/3)U 1 {a)ul(/3)ul(a). This se- 
quence acts over a time 10r to generate the following 
effective Hamiltonian at a single vertex, v. 
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+ X ([<r a^a°} v - 2/3/ 7 [<rW<r*]J + O (0 6 ) 

where [0]„/ p denotes the application of the (up to) 
four-body operator O to the spins meeting at a ver- 
tex v or surrounding a plaquette p; we choose \a\ = 
|0| = | 7 | = 0; and J e = X (l-3^ 2 )/a 7 + 0{4> & ) 
with x = a 2 /3 7 2 (2ft/5r). By repeating the operator se- 
quence a second time with sign reversals a — > — a and 
7 — > —7, we cancel the fourth order terms in 0, giving 
Ui23(~ a., (3, — 7)1/123(0:, /?, 7) that acts for a time 20r to 



generate the effective Hamiltonian 

H£" = J e II °i + X [<7<WV°] v + O (0 6 ) . (3) 

The sequence U^zi.— a, /?, — 7)1/123 (a, /3, 7) is specifically 
designed to cancel the lowest-order (0 4 ) perturbation 
terms without affecting the gap. The remaining c/> 5 term 
is a 2-body perturbation to Hq C . Repeating this se- 
quence with appropriate sign reversals will cancel these 
higher order terms. However, the ground state sub- 
space of Hq is robust to these remaining perturbations 
(see below). A shorter operator sequence may then be 
preferable to reduce gate errors. The plaquette operator, 
fl-xxxx^ can jj e generated by cyclic permutation of the 
Pauli operators in the above expressions for Ei,E2 and 

Simulation of Hj c then requires application of the 
pulse sequence to all vertices and plaquettes conforming 
to a two dimensional square lattice with periodic bound- 
ary conditions. Vertex and plaquette terms may be ap- 
plied serially as: exp(—iH x g xx t/H)exp(—iH z g zz t/h) « 
exp —i (H x g xx + H z § zz ) t/h. Because only the perturba- 
tion terms fail to commute, the truncation error in the 
above expression occurs at orders larger than cf> 7 . For 18 
pseudospins, representing a 3 x 3 system with toroidal 
boundary conditions, a completely serial implementation 
yields a stroboscopic cycle time of 720^s using the es- 
timate t ~ 500ns [l(| and the minimal count of one 
CPHASE and four 1-spin gates to realize all Uj(a) [rH |. 
This may be reduced by implementing some operators in 
parallel. 

Simulated Thermalization .- The pscudospin subspacc 
of the system defined by the internal states of the 
trapped atoms will interact with the external environ- 
ment through the controlled quantum operations in the 
above pulse sequences and uncontrolled noisy interac- 
tions. Noise in the optical lattice system will not drive the 
simulation subspace to a state that is thermal under the 
simulated Hamiltonian 18|. Additionally, noise in the 



above sequence of control gates will add entropy and ef- 
fectively heat the system. The entropy production AS re- 
sulting from imperfect gate operation AS ~ EPG, where 
EPG is the error per gate [14| . Quantum circuit models 
are usually supplemented by error correction schemes to 
effectively remove entropy from the system. We take a 
different approach here, constructing an effective system- 
reservoir interaction to control system entropy and relax 
the system to the ground state or a thermal state. 

To maintain the simulated system at a thermal steady 
state we add an interaction H sr of the system pseudospins 
with a set of ancillary pseudospins. In the optical lattice 
system, these ancillary pseudospins, which may be a 2 nd 
species of atom, are trapped in an offset, intercalated op- 
tical lattice, such that each ancillary atom is adjacent 
to a system atom. Consider a Hamiltonian with local 
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n-body interactions of the form: H = — Yl v Jv StV % 
, where h^f is an n-body operator involving a neighbor- 
hood of pseudospins, A/", including pseudospin i, with 
eigenvalues ±1, v labels the type of interaction, and J v 
is a constant. Additionally we define the pseudospin flip 
operator such that ?%\h v N = ±1) = \h" N = Tl) when 
i £ N. When all [h^j-, hfr,] = 0, as is the case for Eq. <JTJ) , 
we can define E\ v = |E? (1 + hfr) (1 + hfc), T iyV = 
(1 - hf/) (1 + h u H ,) with i = N n TV'. These areare 
(2n — l)-body interactions; E\ v creates a pair of exci- 
tations about i and Ti „ translates an excitation about 
i. The energy gap for creation of a pair of excitations is 
A„ = 4J„. 

A route to guaranteeing the thermal equilibration of 
this system is for it to evolve under the Lindblad mas- 
ter equation p = — i/h[Ho,p] + L[p], where p is the 
density matrix and L [p] is the superoperator L [p] = 
J2u( 2c uP c l- 4 c wP- p4 c o>)> with {c w ,4} the Lind- 
blad operators. With given by 
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, (4) 



the Lindblad master equation describes equilibration 
with a bath of temperature T = — A/ In (p). The unique 
stationary state of the system is then the thermal state 
under Hq with temperature T. A* and 7* are relax- 
ation rates, and their values dictate the thermalization 
time. For simplicity wc have set A v — A. To gen- 
erate evolution under such a master equation, we in- 
troduce a set of non-interacting ancillary pseudospins 
that independently undergo strong dissipation. Each lo- 
cal neighborhood of the system interacts locally with 
a single thermal ancillary pseudospin Tj it , and a single 



maximally mixed ancillary pseudospin 



for each 



Hsr = .'/V„ (/•.; , • i ll o u • //.c.) The master 
equation of the system and ancilla pseudospins combined 
is of the above Lindblad form with Lindblad operators: 



{<H,v} = 



where A and 7 define the relaxation rates of the individual 
ancillary pseudospins. With this choice of H sr it can be 
shown that for g <C A the system pseudospins evolve 
under a rcnormalized Lindblad master equation with c u 
given by Eq. (j4|) and A* = 4(g/ti) 2 A, thus leaving the 
thermal state of Hq as the unique stationary state. 

For p = 0, the effective system- reservoir interaction 
cools the system towards the ground state, and the 
Lindblad operators can be reduced to the n-body terms 
{VX*(Ei, u +T itV ),V¥(E itV +Tl v )} dlEi In this limit 
the ancillary pseudospins become an effective low tem- 
perature bath with a cooling rate T c ~ A* and heat- 
ing rate determined by gate errors and any environmen- 



tal noise. Competition between these rates leads to a 
minimum reachable temperature for the system T m ; n ~ 
A/ln(r c /r c ), where L c - EPG x fi, with EPG and 
the error rate and frequency of application of Uj(a), re- 
spectively. 

The Lindblad master equation, with operators given 
by Eq. ([5]), generates a unitary system-reservior interac- 
tion but nonunitary reservoir relaxation. Stroboscopic 
simulation of H sr is performed in a manner analogous to 
the Hq simulation described above. Phase angles are 
chosen in the 1- and 2-body gates to generate an effective 
static interaction strength g over the time t sr between ap- 
plications of H sr , such that gt sr /h < tt/2. Nonunitary 
evolution of the reservoir is generated by encoding the 
reservoir as two levels of a A-system. The pseudospin 
states are the ground state |0) and the meta-stable state 
|1). State |2) is chosen to have fast spontaneous emission 
to |0), with rate T20. This spontaneous emission is the 
decoherence mechanism required to generate the nonuni- 
tary Lindblad evolution. The ancillary pseudospin levels 
can be placed in a thermal state via the following proce- 
dure: i) 7r-pulse on the |1) — > |2) transition, ii) Wait for 
decay to ground state, |0). iii) 7r-pulse on |0) — > |1) tran- 
sition, iv) 0-pulse on the |1) — > |2) transition, v) Wait 
for decay, which now yields the final pseudospin state, 
p = diag {sin 2 (6*), cos 2 (6*)} , corresponding to an effective 
temperature T e g = A/ (2 In (cot 0)). The above strobo- 
scopic procedure generates A* ~ g 2 t sr /h 2 in Eq. j4|. The 
procedure can be simplified in the limit of cooling towards 
zero temperature by eliminating steps iii-v, when it be- 
comes similar to the optical pumping scheme employed 



in measurement of qubit states for trapped ions 21 1. 
This thermalization procedure is then repeated and in- 
terleaved with the stroboscopic application of Hq . 

Thermalization of the Toric Code.- Hj 1 ^ is in the lo- 
cal form of Hq, with two types of excitations, electric 
charges and magnetic vortices, {y = e, m) that reside on 
vertices and plaquettes, respectively, of the square lattice. 
The excitation operators are defined with h% = Y\j ev c|, 
hp 1 = YljfzpO-j, and E^' m = af' z . Each link must inter- 
act with four ancillary pseudospins in the limit T <C x 
or X "C g to allow thermalization to the ground state or 
the thermal state of Hq C , respectively. 

The stroboscopic generation of Hq C outlined above 
introduces truncation perturbations in the perturbativc 
expansion, e.g, the second term in Eq. ([3]), which are 
distinct from extrinsic errors due to experimental noise 
and gate inaccuracies. If sufficiently large, such trunca- 
tion perturbations could drive the system away from the 
desired ground state phase . We now show on a finite 
sized system accessible to current experiments [13], that 
the intrinsic perturbations can be kept sufficiently small. 



Fig. Ufa) plots the gap of H^ 



HTC 



xJ2(i j) a i aV i as a function of the strength of the pertur- 
bation for a 3x3 planar lattice with toroidal boundary 
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Figure 2: (a) Energy spectrum (in units of J e = Jm) and 
(b) ground state fidelity vs. perturbation strength from exact 
diagonalization of the 18 site toric code with h z = 0.05. 



conditions (Fig. [1} . The Zeeman field is added here to 
fully split the ground state degeneracy and ensure robust 
characterization of the eigenstates of H even in the pres- 
ence of small additional perturbations. We define the 
ground state fidelity as F^ s = |(*° |* n (x, h z ))\, where 
the |^°) are the degenerate ground states of Hj c and 
{Xi h z )) are the nearly degenerate ground states of 
h ■ Fig- Hit") shows the ground state fidelity as a 
function of \- This fidelity determines the robustness of 
topological operations that will be performed via string 
operators to measure or perform gates on the system. 
We see that for \x\ < 0.4 the features of the topological 
phase persist, including the approximate four-fold degen- 
eracy of the ground state and a finite gap to excitations. 
This corresponds to a maximum value of </> ~ 0.4, which 
constrains the gate operations in the pulse sequences, 
Eq. 12]). This robustness should increase with increas- 
ing lattice size, and is consistent with known stability of 
Hq C to perturbations p2| . 

Increasing <p increases J and therefore the gap of the 
Hq C ; however it also increases x/J> which reduce the gap 
of -£T^Y 2 and topological protection for large \- We also 
note that J ~ 1/Nq where Nq is the number of sequen- 
tial gates used to simulate H . For larger lattices, some 
degree of parallelization is thus desirable to ensure that 
the gap does not decrease with the lattice size. Choosing 
X = 0.2, the gap achieved by a completely serial imple- 
mentation is A w 0.6 fiK/N sys , where N sys is the number 
of system atoms used. With the cooling sequence seri- 



ally interleaved^ « 0.1 ^K/N syB and A* - 10 4 s~ 1 /N sys 
are achievable jl4j . For a minimal system of 18 system 
atoms, this allows for an effective temperature T c ff < A 
to be reached with an error rate of EPG ~ 1CD 4 or 
less [14 1- 

Sources of Errors.- This scheme is designed to be ro- 
bust against errors within the pseudospin subspace. The 
dominant source of residual error in the implementation 
discussed here is leakage from the Rydberg levels due to 
spontaneous emission and black body radiation. The lat- 
ter may be effectively suppressed by working at low tem- 



peratures [23j | . and spontaneous emission is minimized 
by utilizing states with larger n. With n < 180, we es- 
timate that spointaneous emission errors can be reduced 
to ^10 -6 per gate, allowing for up to 10 3 stroboscopic 
cycles. 

Discussion - We have developed a formalism for the stro- 
boscopic generation of n-body Hamiltonians using 1- and 
2-body quantum operations together with a dissipativc 
thermalization scheme. We have applied this here to the 
toric code Hamiltonian in the context of addressable op- 
tical lattice experiments [13]. Our method applies to 
a wide range of lattice spin models [3] as well as to 
other experimental setups 24J. The dynamic generation 
both of a Hamiltonian possessing a topologically ordered 
ground state and of an effective thermalization mech- 
anism offers the possibility of robust simulation of the 
ground state and of the creation and braiding of anyonic 
excitations [25j . These are essential components required 
for the topologically protected storage and manipulation 
of quantum information. 

We thank D. Weiss for useful discussions. This ma- 
terial is based upon work supported by DARPA un- 
der Award No. 3854-UCB-AFOSR-0041. During the 
preparation of this manuscript related results discussing 
ground state preparation were reported 26 1. 
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